Critical properties of a continuous family of XY noncollinear magnets 
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Monte Carlo methods are used to study a family of three dimensional XY frustrated models 
interpolating continuously between the stacked triangular antiferromagnets and a variant of this 
model for which a local rigidity constraint is imposed. Our study leads us to conclude that generically 
weak first order behavior occurs in this family of models in agreement with a recent nonperturbative 
renormalization group description of frustrated magnets. 

PACS numbers: 75.10.Hk, 75.40.Cx, 75.40.Mg 



In spite of intensive study during the last 25 years the 
critical behavior of XY or Hcisenberg frustrated magnets 
with a noncollinear ground state is still a strongly de- 
bated topic (see Refs. P, @ and references therein) . This 
is, for instance, the case of the celebrated XY Stacked Tri- 
angular Antifcrromagnct (STA) that we consider here, 
whose Hamiltonian is: 

" »..s,.s ; (1) 
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where the Si are two-component vectors and the sum 
runs over all pairs of nearest neighbor spins on a stacked 
triangular lattice. In Eq. the in-plane interactions 
are antiferromagnetic (J < 0) and the inter-plane in- 
teractions are taken to be ferromagnetic (J > 0) with 
| J | = 1. The competition between the in-plane antiferro- 
magnetic interactions produces a ground state where the 
three spins Si, S2 and S3 on each elementary triangular 
plaquette are oriented at 120° one to another, that is: 

S! + S 2 + S 3 = 0. (2) 

For such a system, the order parameter is given by two 
orthogonal vectors of the same norm, a fact that has 
led to the hypothesis that noncollinear magnets could 
undergo a second order phase transition associated with 
a new — chiral — universality class Q . 

However, it is now becoming more and more widely 
accepted that the physics of frustrated magnets is not 
given by such a simple picture 00, IE 01 The complexity 
of the critical behavior of STA systems is exemplified in 
the experimental results (for a review see Ref. and 
references therein). Indeed, scaling behavior is generally 



found while universality is violated since materials which 
belong, a priori, to the same universality class display 
different critical exponents. It is thus difficult to interpret 
the experimental data within the usual picture of critical 
phenomena. 

On the theoretical side, the situation is controver- 
sial 0, 0- A nonperturbative renormalization group 
(NPRG) method 0, 0] finds that the phase transitions 
in these systems are generically weakly first order with 
the possibility of standard strongly first order transitions. 
More precisely while there is no RG fixed point there ex- 
ists a whole region where the RG flow is very slow, cor- 
responding to large — but finite — correlation lengths 
at the transition, which allows one to compute pseudo- 
critical exponents. Hence the violation of universality 
observed experimentally has a natural explanation since, 
in the absence of a fixed point, the pseudo-critical ex- 
ponents found depend on the microscopic Hamiltonians. 
On the other hand, Pelissetto et al. have derived and re- 
summed the six- loop (3 functions in three dimensions • 
They find a fixed point and therefore predict a second 
order phase transition. Calabrese et al. have claimed 
that the lack of universality found can be accounted for 
by the spiral nature of the approach to the focus fixed 
point 0- 

Our aim is to discriminate by means of numerical sim- 
ulations between these two scenarios that both predict 
scaling behavior with varying critical exponents but dif- 
fer in the predicted asymptotic behavior for tempera- 
tures very close to the transition temperature. Different 
systems having the same symmetry breaking scheme as 
XY-STA but differing by microscopic details have been 
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studied numerically 

S H El El 111 13 G3- Apart 
from STA, they have all been found to undergo first order 
transitions @,HE1E|] ^ s ^ or STA itself earl y nu merical 
simulations favor a continuous transition |l2j, |l3j, LLJ, LL5j • 

However the most recent simulation performed with lat- 
tice sizes much larger than those considered previously 
leads to first order behavior Thus, although there 

are now strong indications in favor of the nonperturba- 
tive scenario, the numerical evidence is still too weak to 
lead to a definite conclusion. In order to shed some light 
on this problem we consider a family of models, built up 
as generalizations of STA. We find that they all exhibit 
cither strong first order transitions or a scaling behavior 
with varying exponents associated with a weak first order 
phase transition. 

To perform our study we partition the original trian- 
gular lattice — with lattice spacing a — into elementary 
plaquettes labelled by an index /. They are composed 
of three spins S{, S2 and S3. The super lattice made of 
these plaquettes is, again, a triangular lattice, with lat- 
tice spacing \fia. The family of models we consider is 
defined by: 

H ( r ) = - E J y S ^ + r E + S 2 + S a) 2 (3) 

m 1 

where the first term is the same as in Eq. Q and where 
the sum over / in the second r-term runs over the indices 
of the super lattice. For r < — 0.5|J|, the ground state 
has the spins aligned ferromagnetically within a plaquc- 
tte but are oriented at 120° from plaquette to plaquette; 
this is not the ground state expected for real materi- 
als. For r = — 0.5| J\ the interactions inside the plaque- 
ttes vanish and the lattice becomes effectively a stacked 
Kagome lattice. For r > — 0.5|J| the ground state and 
thus the symmetry breaking scheme and order parameter 
are identical to those of the STA model, Eq.QJ; we thus 
focus on this range of values in the following. For any 
value of r > the r-tcrm favors locally, that is inside 
the plaquettes, the 120° ground state configuration. It 
thus penalizes the relative fluctuations of the spins S{, Sj 
and S3 inside a plaquette but still lets the orientation 
of the three-spin structure freely fluctuate from plaque- 
tte to plaquette. Note that, since the fluctuations of the 
spins inside a plaquette are noncritical, the r-term should 
not alter the critical physics of the STA if it is universal 
- see below. The r — -> 00 limit consists of a STA for 
which the ground state 120° structure is locally imposed 
at all temperatures. This model is called the STAR (R 
for rigid). Continuous changes in the 'rigidity parameter' 
r from zero to infinity thus correspond to a continuous 
change from the STA to the STAR. The STAR has been 
studied by Loison and Schotte ^3 an d nas been found 
to undergo a strong first order transition. This already 
shows that, in fact, the critical behaviors of STA- like sys- 
tems strongly depend on the microscopic details of the 
model. 
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TABLE I: Critical exponents for various r. The results for 
r = were obtained in [l3|l . 



We use a classical Monte Carlo Metropolis algo- 
rithm [lij to study this generalized model as a function of 
temperature and the rigidity parameter r. A number of 
different thermodynamic quantities have been computed 
with < r < 8 for lattices of linear sizes L = 18 to 
L = 138. In order to reduce the critical slowing down 
at the transition, each simulation run starts from a dis- 
ordered state at a temperature well above the transition 
temperature and is slowly cooled to low temperatures. 
The system is allowed to equilibrate for 3 x 10 5 Monte 
Carlo sweeps (MCS) and the averages were performed 
over 5.5 x 10 5 steps. In order to avoid correlations be- 
tween subsequent measurements, 20 — 30 Monte Carlo 
sweeps were performed between measurements. 

< r < 1.0: the weak first order region. We begin 
our discussion with the small r region which is the most 
interesting one since real compounds arc supposed to be 
close to STA. Moreover, for some of these compounds as 
well as for STA with r = scaling behavior is found 
Since scaling laws are found to hold in this range of r, 
we analyze them by means of the standard methodology 
used for second order phase transitions. The fourth or- 
der cumulant Um of the order parameter displays behav- 
ior characteristic of a second order transition and can be 
used to estimate the transition temperature T*. We have 
determined T* , the order parameter critical exponent /?, 
the correlation length exponent v, the susceptibility ex- 
ponent 7 and the anomalous dimension exponent 77 using 
finite size scaling (FSS) methods [l3j]. For appropriate 
values of T* and of (3 and v, the data should collapse 
onto a single curve. A sample log-log finite size scaling 
plot of the order parameter M for r = 1.0 is shown in 
figure^ A similar FSS analysis of the susceptibility also 
satisfies scaling and can be used to determine 7. For 
sizes L = 48, 60 there are some deviations from scaling 
at small values of the reduced temperature t. 

At T* , the structure factor depends on system size as 
S(Q, L) ~ L 2 ~ v where the exponent 77 is the anomalous 
dimension and Q is the ordering wave vector for the three 
sublattice structure. Throughout the entire range < 
r < 1.0, we again find that scaling is satisfied for system 
sizes L < 42. Tabled summarizes our results for the 
critical exponents in the range < r < 1.0. The critical 
exponents for which finite size scaling holds change with r 
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FIG. 1: Log-log finite size scaling plot of order parameter 
M for the r — 1.0 system. The transition temperature is 
T* = 1.7369 and the exponent f3 = 0.2134. The best collapse 
of data was obtained for v — 0.4817 . 



FIG. 2: The dashed line corresponds to a linear fit of the 
energy difference of peak maxima for r = 1,2,4 and 6. The 
point at r = was estimated from the histograms in figure 3. 



showing clearly that universality is violated. Moreover r\ 
is always found to be negative — sometimes significantly, 
e.g. for r = 0.2 — which is impossible if the transition is 
of second order These results strongly suggest that 
scaling aborts ultimately — that is, very close to T* - 
for all r € [0, 1] and that, therefore, the transitions are 
all of first order. 

r > 1.0: the strong first order region. We now con- 
sider the region r > 1.0. As expected, the first order 
behavior found for the STAR, that is for r — ► oo, per- 
sists for finite r. We have been able to follow it down 
to r ~ 1.0 by studying the probability distribution for 
the energy P(E) . This probability distribution is a use- 
ful quantity to locate first order transitions. Away from 
the transition point, one expects a Gaussian distribution 
in energy while at the transition point a double peak in 
P(E) occurs if it is first order. This point can be located 
by finding the temperature where the weights under the 
peaks in P(E) become equal. However distinguishing be- 
tween a weak first order and a second order transition is 
rather delicate especially in case of a possible tricritical 
point. In order to make sure what type of phase transi- 
tion the double peak structure represents, it is necessary 
to study the size dependence of the probability distribu- 
tion P(E). In the case of a first order transition, the 
separation of the peak energies E\ and Ei equals the la- 
tent heat up to 1/L corrections Thus two peaks 
remain separated with increased size L and, as shown in 
[lij . become sharper. 

We have found a double peaked structure in P(E) for 
rigidity parameters r = 1.0,1.5,2.0,4.0,6.0 and 8.0 for 
linear sizes L < 60. As the value of r decreases, the size 
of the system where the double peak structure starts to 
appear increases. However, the energy autocorrelation 
time is less than 2000 MCS in this range. The dou- 



ble peaked structure in P(E), and thus the first order 
behavior, is apparent but only for lattice sizes L > 48 
for r = 1.0. If we extrapolate to the thermodynamic 
limit, we can estimate the latent heat AE = \E\ — E%\ 
at the transition. Using the values of AE for r > 1.0 
and extrapolating to the more interesting small r region 
predicts that AE remains nonzero at r = with an es- 
timated value of 0.010 ± 0.002. This behavior suggests 
that the phase transition at r = might be weakly first 
order. 

In figure El we plot AE vs r where it can be seen that 
it extrapolates linearly to the physically interesting limit 
of small r. The maximum size used for r = 6 was L = 36 
and hence the errorbars are larger than those at smaller 
r. The value at r = is finite indicating a weak first or- 
der transition. The extrapolation also predicts that AE 
vanishes at ro = — 0.47± 0.02. Our error bars for ro are 
compatible with r = —0.5, corresponding to the stacked 
Kagome lattice, for which one expects a completely dif- 
ferent physics. We conclude that STA undergoes a very 
weak first order phase transition and that there is no 
tricritical point for all systems described by H(r) with 
r > -0.5. 

STA r = 0. In order to confirm our result of a non- 
zero latent heat at r = 0, we have studied the energy 
probability distribution for the r = case using much 
larger system sizes. Figure 01 shows P{E) at T* for sizes 
L = 87, 111 and 138. A double peak structure charac- 
teristic of a first order phase transition is now evident 
at these sizes. As the size of the system is increased, 
the peaks sharpen and the depth of the minimum be- 
tween the peaks also increases. For these larger sizes the 
energy autocorrelation time increases significantly to val- 
ues greater than 12000 MCS. This critical slowing down 
could perhaps be improved using ovcrrclaxation or multi- 
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canonical histogram methods. [18( However, we estimate 
a value AE ~ 0.013 ± .001 from figure |3] which is in excel- 
lent agreement with the previous estimate 0.010 ± 0.002 
obtained from the extrapolation of AE(r) from larger 
r > 1.0 to r = 0. 

In conclusion, we have studied a system interpolating 
continuously between the XY STA and STAR models 
and find unambiguous evidence that the phase transi- 
tions in these systems is either of first order or weakly 
first order accompanied by scaling with exponents vary- 
ing with r. These results provide strong arguments in 
favor of the nonperturbativc approach to the physics of 
frustrated magnets Q and suggest that, even if a fixed 
point exists, its relevance for STA-likc compounds whose 
microscopic Hamiltonians are supposed to be well de- 
scribed by Eq. Q), is doubtful. Indeed, while the two 
— perturbative and nonperturbative — approaches both 
predict scaling in a transient regime near T* , they differ 
in several ways. First, asymptotically, i.e. at T* — the 
scaling behavior aborts in the nonperturbative scenario 
while it persists in the perturbative one since it is associ- 
ated with a true continuous phase transition. Second, the 
weak first order transitions are generic in the nonpertur- 
bative scenario whereas they should be exceptional in the 
perturbative one Q. Indeed, in the presence of a fixed 
point, systems undergoing weak first order phase tran- 
sitions correspond, in the space of coupling constants, 
to initial conditions of the RG flow that are outside the 
basin of attraction of the fixed point but that lie in the 
vicinity of the separatrix between the first and second or- 
der regions. This behavior should not be generic since it 
implies a fine-tuning of the parameters of the microscopic 
Hamiltonian. Note however that one cannot completely 
exclude that, by varying another microscopic parameter 
like the ratio between the interlayer and intralaycr cou- 
pling constants, the system be driven into the basin of 



attraction of a fixed point. 

The present work as well as the one performed by 
Itakura [TlT | predict a first order phase transition for 
STA. Previous numerical simulations failed to reach the 
asymptotic regime where the correlation length satu- 
rates, that is, where the scaling aborts. This suggests 
that the same phenomenon could occur in experiments, 
this time due to the difficulty to probe the close vicinity 
of T* , and that all real materials also display weak first 
order transitions. This hypothesis is reinforced by the 
fact that the exponents found in our study are close to 
the experimental ones for STA compounds (see and 
references therein). Let us consider for instance the ex- 
ponent (3 which is the best measured one. Here we find 
a variation of (3 in the range [0.213,0.253] while for XY 
STA compounds one has, for instance, (3 ~ 0.228 for 
CsMnBr 3 , (3 ~ 0.243 for CsNiCl 3 and f3 ~ 0.23 - 0.25 for 
CsCuCbj which has been finally found to undergo a first 
order transition when the study of the close vicinity of 
T* has been refined. 

An important open question concerns the Heisenberg 
case where there exists the same conflict between the 
perturbative and nonperturbative scenarios. However, 
according to the NPRG scenario, the transition should 
be even more weakly of first order for Heisenberg than 
for XY spins so that the controversy will be more difficult 
to settle. 
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FIG. 3: P(E) for STA XY model (r = 0) for L = 87, 111 and 
L = 138 at T*. 
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